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Abstract 

We prove the global strong solvability of a quasilinear initial-boundary value problem 
with fractional time derivative of order less than one. Such problems arise in mathematical 
physics in the context of anomalous diffusion and the modelling of dynamic processes in 
materials with memory. The proof relies heavily on a regularity result about the interior 
Holder continuity of weak solutions to time fractional diffusion equations, which has been 
proved recently by the author. We further establish a basic L2 decay estimate for the special 
case with vanishing external source term and homogeneous Dirichlet boundary condition. 
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1 Introduction and main results 

Let T > 0, N > 2, and SI C be a bounded domain. The main purpose of this paper is to 
prove the global strong solvability of the time fractional quasilinear problem 

df{u - u ) - D^i^Dju) = f, t G (0,T), x € fi, 

u = g, te(0,T),xeT, (1) 
u\ t=Q = u , x G fi, 

where we use the sum convention. Here Y = dQ, Du denotes the gradient of u with respect to 
the spatial variables and df stands for the Riemann-Liouvillc fractional derivation operator with 
respect to time of order a G (0, 1); it is defined by 

d"v(t,x) = dt / gi- a (t — t)v(t, x) dr, t > 0, x G O, 
J o 

where gp denotes the Riemann-Liouville kernel 

t fs-i 

The functions / = f(t,x), g = g(t,x), and uq = uq(x) are given data. 
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During the last decade there has been an increasing interest in time fractional diffusion 
equations like ([TJ and special cases of it. An important application is the modelling of anomalous 
diffusion, see e.g. the surveys [TS], [TB]. In this context, equations of the type ((TJ) are termed 
subdiffusion equations as the time order a is less than one. While in normal diffusion (described 
by the heat equation or more general parabolic equations), the mean squared displacement of 
a diffusive particle behaves like const-i for t — >• oo, in the time fractional case this quantity 
grows as const-f", for which there is evidence in a diverse number of systems, see |15j and the 
references therein. In the case of equation ([T]) the diffusion coefficients are allowed to depend on 
the unknown u. 

Another context where equations of the type (JTJ) and variants of them appear is the modelling 
of dynamic processes in materials with memory. An example is given by the theory of heat 
conduction with memory, see [19] and the references therein. Another application is the following 
special case of a model for the diffusion of fluids in porous media with memory, which has been 
introduced in pQ: 

8f{q - q ) - div (K(q)Dq) = /, t G (0, T), x £ Cl, 

q = 0, <e(0,T),xer, (2) 
q\t=o = qo, x e fi. 

Here a £ (0, 1), q = q(t,x) denotes the pressure of the fluid, k = n(q) stands for the permeability 
of the porous medium, and / is related to external sources in the equation of balance of mass. 
Model ^ is obtained by combining the latter equation with a modified version of Darcy's law 
for the mass flux J which reads 

J=-d}- a {K(q)Dq), 

and by assuming that the (average) mass of the fluid is proportional to the pressure. We refer 
to [9] , where a more general model is discussed. 

We next describe our main result concerning ([1]) . Letting p > N + ^ we will assume that 

(Ql) / G L p ([0,T];L p (n)), g G ^"^([O, T]; L p (T)) n L p ([0, T]; B„> (r)), u G S^^(O), 
and u = g\ t=0 on T; 

(Q2) A = (aij)i,j=i,...,iv S C 1 (R; Sym{A^}), and there exists v > such that ^ H£| 2 

for all y e'l and £ G R N . 

Here Sym{A} denotes the space of A^-dimensional real symmetric matrices. For s > and 
I < p < oo the symbols and B^ p refer to Bessel potential (Sobolev spaces for integer s) and 
Sobolev-Slobodeckij spaces, respectively. 

The main result concerning the problem (J1J reads as follows. 

Theorem 1.1 Let fl C R N (N > 2) be a bounded domain with C 2 -smooth boundary T. Let 
a 6 (0, 1), T > be an arbitrary number, p > N + ^, and suppose that the assumptions (Ql) 
and (Q2) are satisfied. Then the problem {Ip possesses a unique strong solution u in the class 

u G H«([0,T};L p (Cl))nL p ([0,T};H%(n)). 

Note that short-time existence of strong or classical solutions to problems like ([lj can be estab- 
lished by means of maximal regularity for the linearized problem and the contraction mapping 
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principle. This has been known before, see e.g. @], ED' an d for the semilinear case [5]. For 
results on maximal regularity for fractional evolution equations we further refer to [3] , [12] , [19] , 
and [35]. The novelty here is that T > can be given arbitrarily large without assuming any 
smallness condition on the data. In some papers (global) generalized solutions are constructed 
for quasilinear subdiffusion problems, see [7] and [TU]. These results are based on the theory of 
accretive operators. 

The crucial step in our proof of the global existence result is an estimate of the Holder norm 
of u on parabolic subdomains [0,6] x Cl which is uniform with respect to S £ (0,T]. In a very 
recent work ([53], see also [21]), the author was able to prove interior Holder regularity for weak 
solutions of time fractional diffusion equations of the form 

d?{u - u ) - D i (a ij (t,x)D j u) = /, t £ (0,T), x £ n, (3) 

with a £ (0, 1) and merely bounded and measurable coefficients, see Theorem 12.21 below. Using 
this result, we derive conditions which are sufficient for Holder continuity up to the parabolic 
boundary. Here we do not aim at high generality but we are content with finding some simple 
conditions which are also necessary when studying the quasilinear problem (fT]) in the setting of 
maximal L p -regularity. Because of the latter it is natural and also not so difficult to use the 
method of maximal L p -regularity to achieve the goal. 

In this paper we also prove a basic decay estimate for the L2(0)-norm of the solution u to 
([TJ in the special case when / = g = 0. It is shown that for the global strong solution of (JTJ we 
have in this case 

KV)li 2( n) < t>0, (4) 

where c = c(a) and /i = n{y, N, fi) are positive constants. A polynomial decay estimate has 
already been known in the linear case for problems of the type 

df{u - m ) - D l (a ij {x)D J u) = 0, t £ (0,T), ie!!, 

n = 0,te (0,T), xeT, (5) 
u\ t= o =u , x e n, 

see [T7] Cor. 4.1]; for the special case a%j = Sij we also refer to [TJ] p. 11-13]. By means of 
an eigenfunction expansion for u, there it is shown that in general \u(t, ')|i 3 (n) decays as l/t a 
as t — > oo, which is optimal w.r.t. the exponent and stronger than ((4]), where we have this 
behaviour for the square of \u(t, ■)\l 2 (Q)- I n particular we do not have exponential decay as in 
the case a = 1. 

The smaller exponent (a/2 instead of a) in the quasilinear case is due to the different method. 
Note that the aforementioned method is no longer applicable in the quasilinear case. Our proof 
of ([4]) is based on energy estimates; it makes use of the fundamental convexity inequality ([9]) for 
the Riemann-Liouville fractional derivative, see Theorem 12.41 below. Our method applies to a 
more general class of quasilinear problems, including e.g. the time fractional p-Laplace equation, 
and it extends to the (more natural) weak setting. This will be elaborated in a forth-coming 
paper. 

The paper is organized as follows. In Section 2 we collect some known results on a priori 
estimates for linear time fractional diffusion equations and recall the fundamental convexity 
identity (j9]) for the fractional derivative. Section 3 is devoted to Holder regularity up to t = 
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for weak solutions of Regularity up to the full parabolic boundary is then established in 
Section 4. Using these estimates, the global existence result, Theorem ll.il is proved in Section 
5. Finally, we derive the described decay estimate in Section 6. 



2 Preliminaries 



We first fix some notation. For T > and a bounded domain £1 C 1^ with boundary L we put 
il T = (0,T) x Q and T T = (0,T) x T. The Lebesgue measure in R N will be denoted by Ajv- 

The boundary L is said to satisfy the property of positive geometric density, if there exist 
/3 G (0, 1) and po > such that for any xq G L, any ball B(xo,p) with p < pg we have that 
\ N (ttn B(x ,p)) < p\ N (B(x ,p)), cf. e.g. Section LI]. 

By y + :— max{j/, 0} we denote the positive part of y G K. 

In the following we collect some known results from the linear theory which are basic to the 
investigation of (P). 

Let T > and il be a bounded domain in with N > 2. We consider the linear time 
fractional diffusion equation 



d?(u - u ) - D l (a lJ (t,x)D J u) = /, t G (0,T), x G ft. 

We assume that 
(HI) Ae L 00 {il T ;R NxN ) 7 and 

N 

\aij(t,x)\ 2 < A 2 , for a.a. (t,x) G fir- 

(H2) There exists v > such that 

aij(t,x)£,i^j > ^|£| 2 , for a.a. (t, x) G ftr, and all £ G 

(H3) u Q G Loo(^); / 6 L r ([0,T];Lg(Q)), where r,?> 1 fulfill 

1 TV 

— + 5- = 1 - 

ar 2q 



(0) 



and 



r G 



1 



a(l — k) 



, oo 



, q e 



JV 



2(1 - k) 



, oo 



K G (0,1). 



Following |23 j and [25] we say that a function u is a weafc solution of (J6j) in f^, if u belongs to 
the space 

S a := {v G L 2 /(i-a),«,([0,T];i 2 (fi)) n L 2 ([0,T]; ff^fi)) such that 
5 i_ Q * w G C([0, T];L 2 (fi)), and ( 5l _ a * u)| t=0 = 0}, 



and for any test function 

^ffj'N :=i/ 2 1 ([0,T];i 2 (r!))nL 2 ([0,T];i7 2 1 (fi)) (^(Q) :=Q^(Q)^ 



(n) 
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with r]\t—T = there holds 

rT 



I I ( - Vt[gi-a *(u- uo)] + aijDjuDirj) dxdt = I / frjdxdt. (7) 
Jo Jo v y Jo Jn 



Here L Pj stands for the weak L p space and f\ * f-2 means the convolution on the positive halflinc 
with respect to time, that is (/i * f2){t) — f fi(t — 7~)/2(t) dr, t > 0. 

Existence of weak solutions of ([5]) in the class <S Q has been shown in [28 . For example, 
assuming (HI), (H2), uo £ Z,2(f2), an d / £ ^(^t), the corresponding Dirichlet problem has 
a unique solution in S a . Global boundedness of weak solutions has been obtained in [24], the 
result can be stated as follows, cf. [24J Corollary 3.1]. 

Theorem 2.1 Let a E (0,1), T > and Q be a bounded domain in R N (N > 2) with T 
satisfying the property of positive density. Let further the assumptions (H1)-(H3) be satisfied. 
Suppose that u £ S a is a weak solution of fSJ) such that \u\ < K a.e. on Tt (in the sense 
that (u — K) + , (—u — K)+ E L2QO, T]; H\ (ft))) for some K > liioU^rn) • Then u is essentially 
bounded in Q,t and 

\u\ Loo{ n T )<C(l + K), 
where the constant C = C(a,r,q,T, N,is,Q,\f\ Lr ( Lq }). 

An interior Holder estimate for bounded weak solutions of © has been proved recently in [23] , 
see also [25]. For Pi, /3 2 E (0, 1) and Q C fir we put 

r 1 r \ u {t> x ) ~ u ( s j y)\ \ 

MC31 " 2(Q) " (t,*UsJZ^(*,v) \\t-8\K + \x-v\KJ- 
Then the interior regularity result reads as follows, cf. [23] Theorem 1.1]. 

Theorem 2.2 Let a E (0,1), T > and Vl be a bounded domain in R N (N > 2). Let the 
assumptions (Hl)-(HS) be satisfied and suppose that u £ S a is a bounded weak solution of {TjJ) in 
Q.T- Then there holds for any Q C fir separated from Tt by a positive distance d, 



C'(|u|L oc (n T ) + \u \ Loo (n) + \f\L r ([0,T];L q (Q))J 



with positive constants e = e(A, v, a, r, q, N, diamQ, mii T Z \^Q r) and C = C(A,v,a,r,q,N, 
diamCl, \N+i(Q),d). 



Remark 2.1 The statement of Theorem 12.21 can be extended to the case where the right-hand 
side of equation has the form 

kf kg 
k=l k=l 

with h E L rk ([0,T];L qh (Q)), k = l,...,k f , E^M? £ L rW ([0, T]; L qW (O)), k = l,...,k g , 
and all pairs of exponents (rk,qk) and (r^ k \q^), respectively, are subject to the condition in 
(H3). This follows from [23] Remark 6.1], see also [25]. 
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Next we are concerned with maximal L p -regularity for the corresponding problem in non- 
divergence form, 

d?(u - u ) - a lJ (t,x)D i D j u = f, t e (0,T), x £ ft, 

u = g, tG(0,T),xGT, (8) 
u\ t =o =u , x e Cl, 

where we again use the sum convention. The following result is a special case of |27[ Theorem 
3.4] on linear initial-boundary value problems in the context of parabolic Volterra equations. 

Theorem 2.3 Let T > and £1 C M. N be a bounded domain with C 2 -boundary T, and N > 2. 
Let a € (0,1) and p > — + ^. Suppose that A — (a%j)i,j=i,...,N & C([0,T] x r2;Sym{7V}), and 
there exists v > such that aij(t,x)£i£j > i^l 2 /or aZZ (t,x) £ [0,T] x fi and £ £ R^. Tften the 
problem has a unique solution u in the class 

Z := ff p Q ([0,T];£ p (fi)) nL p ([0,T];ff p 2 (l])) ^ C([0,T] x fi) 

if cmd on/?/ «/ t/ie following conditions are satisfied. 

(i) f e L p ([0,T];£ p (n)), 3 £ F D := ^"^([O, T]; L P (T)) D L p ([0, T]; fl£r * (r)), and u £ 

fijj Uo = g| t=0 n T. 
In this case one has an estimate of the form 

\u\z < C(|/|z, p (n T ) + \g\v D + K|y 7 ), 
where C only depends on ot,p, N, T, CI, A. 

We conclude this preliminary part with an important convexity inequality for the Riemann- 
Liouville fractional derivation operator, which will be needed in Section [6] 

Theorem 2.4 Let a £ (0, 1), T > and % be a real Hilbert space with inner product {-\-)h- 
Suppose that v £ L2([0,T];%) and that there exists x £ T-L such that v — x £ qH^HO, T]; H) := 
{g a * w : w £ £a([0, T]; %)}. TTien 

(«(*), ^ (ffi-a * «)(t)) w > \j f (gi- a * + \ 9l - a (t)\v(t)\ 2 n , a.a. t £ (0, T). (9) 

Proof. This follows from Theorem 2.1, Proposition 2.1, and Example 2.1 in [32]. □ 

3 Regularity up to t = 

The objective of this and the following section is to find conditions on the data which ensure 
Holder continuity up to the parabolic boundary for weak solutions of the linear time fractional 
diffusion equation (J6j) . As already mentioned in the introduction we do not aim at great generality 
but at results which are sufficient for the quasilinear problem to be studied. 
We first discuss regularity up to t = 0. 



G 



d?(w - u ) ™ Aw = 0, te (0, 1), x G 1^, 

JV 



Theorem 3.1 Let a G (0,1), T > and f2 6e a bounded domain in R. (W > 2j. Lei i/ie 
assumptions (Hl)-(HS) be satisfied. Let further Q' d Q be an arbitrary subdomain and assume 
that 

2—2- ~ 1 TV 

uo|q 6 B pp pa (Q) with p>— + —, 

for some C 2 -smooth domain fl such that fi' C Q C f2 and 11' is separated from dfl by a positive 
distance d. Then, for any bounded weak solution u of 0) in Qt, there holds 

H C f:'([0,T]x!F) ^ C(\u\ Laa{nT) + luoU.Cn) + |«o| a _A + \f\L r (lO,T]:L q (Cl))) (10) 

with positive constants e — e(A, z/, a, r, iV, diamft) and C — C(A, ^, Q;,p, r, iV, diam^l, T, 

Ajv(n'))- 

Proof. The basic idea of the proof is to extend u to [— 1 , T] x f2 such that u is Holder continuous 
on [—1,0] x Q' and to apply Theorem 12.21 

To this purpose we first extend uq\q G _B PP p ° (il) to a function zio G i?pp P ° (1^). By [27] 
Theorem 3.1], the problem 

u>|t=0 = uo, a; G 
possesses a unique solution w in the class 

Z := H«(10, 1};L P (R N )) n L„([0, 1]; 13* (R*)), 
and one has an estimate of the form 

Ms ^ Co|"o| »— a. < C \u \ a— a. . • 

B„,*>«(R*) B pp pa (U) 

Note that by the mixed derivative theorem (cf. [21]). 

Z^ J ff^ 1 - ; )([0,l]; J ff^ (M iv ))) CG [o,l], 

and thus Z ^ BUC 5 ([0, l]xR N ) for some sufficiently small 6 £ (0, a/2). In fact, the assumption 
p > i + -j ensures existence of some <; G (0, 1) with a(l — <;) — | > (5 and 2$ — y > J. 

Multiplying w by a suitable smooth cut-off function tp(t) we can construct a function w G Z 
with u>|t = o = uo and w)|t=i = 0. We then extend u to [— 1,T] x by setting u(t,x) = w(—t,x) 
for t G [-1,0) and i£0. 

Next, we shift the time by setting r = t + 1. Put u(r, x) = u(t — 1, x), r G (0, T + 1), x G f2. 
Define further 

,g:=^u-Au, t G (0, 1), x G ft. 

Then g G L p ([0,l] x Q), since u| T6(0 ,i) G H p Q ([0, 1]; L p (Q)) n L p ([0, 1]; H p 2 (ft)) and n| T=0 = 0. 
Furthermore we have for any test function r\ G iJ^QO, T + 1] x f2), 



/ / ( — ?/ T (gi-Q. * u) + DjuDiT] ) dxdr = gr\ dx dr — r)(gi- a * u) dx 

Jo Jq ^ ' Jo Jo. Jh 



(11) 
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On the other hand, we have for a. a. (t, x) G (1, T + 1) x Cl, 

{gi- a *u)(t,x) = (gi- a *u)(t- l,x) + / gi- a (r - a)u(a, x) da 

Jo 

= (.gi-a *(U- U )) (r - 1, x) + 92-a(T)u {x) 

+ gi- a (r - a)(u(a,x) -u {x)) da. 
Jo 

Set 

h(r,x) = gi- a (r)u (x) + / gi- a (r - a)(u(a,x) - u Q (x)) da =: hi(r,x) + h 2 {r,x), 



hij (r, x) = aij (t — l,x), and f(r, x) = /(r — 1, x) for (t, a;) G (1, T + 1) x Cl. Since w, is a weak 
solution of © in J7t, we thus obtain after a short computation that for any 77 G if 2 ' ([0, T+l] xO) 
with ?7| r= T+i = 



T+l 



~'7t * u) + aijDjuDi-qj dx dr = 

/ (f + h)r]dxdT+ / r](gi~ a *u)dx . (12) 

7i in in T=1 



Adding (fTTj) and (JT2J) shows that u is a weak solution of 

d°u - DiikjDju) = /, t e (0,T + 1), x e f2, 

where 

bij (r, x) = X[o,i] (r) + X(i,t+i] (r)ay (r, x) 

and 

/(t, a;) = X[o,i] (r)g(r, x) + X(i,T+i] (t) (/ + (t, x). 

Evidently, X[o,i]( T )s € L p ([0,T+l]xfl) and X(i,t+i](t)/ 6 L r ([0,T+l];L g (fi)). Concerning 
the h-term we clearly have X(i,T+i\{ T )h\ G Loo([0,T+ 1] x SI). To estimate X(i.t+i] { T )h2, we 
employ the Holder estimate 

\u(a, x) - u (x)\ = \u(a, x) - u(l, x)\ < Ci(l - a) s , a £ [0, 1], x G Cl, 

which results from the embedding Z <—> BUC s ([0, 1] x R N ) and the construction of it. It follows 
that forl<r = t + l<l + T and x G Cl 



|Aa(T,a:)| < Ci / [-ffi-a(r - a)](l - a) s da 
Jo 

[ (t + a)- 1 -'^ 6 da. 
Jo 



r(l-a) 
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Assuming that t = r — 1 G (0, 1) we then have 

Mt,x)\< ( [\t + a)- 1 - a a s da+ [ (t + a)' 1 -** a 5 da 

r(i - a) V7 7 t 



< ad ^ a+s fl 



r(l — a) \a a — S > 

This shows that X(i,t+i] ( T )^2 S L ro ([0,T + l];Xoo(f2)) for all 1 < r < — In particular we 
find some f > — such that X(i,t+i]( t )^2 S £f ([0, T + 1]; ioo(fi)). 

All in all we see that / is of the form / = Ylt=i /«' wnere fi € in([0, T + 1]; (fi)) with 

1 N 
+ — <1, i=l,2,3,4. 

Hence Theorem l2.2l and Remark |2 . 1 1 imply that u is Holder continuous in [1/2, T + 1] x fl'. This 
in turn yields Holder continuity of u in [0,T] x f2', and it is not difficult to see that u is subject 
to the estimate (1101). □ 



Remark 3.1 It follows from Remark 1 2. II and the proof above, that Theorem 13 . 1 1 can be gener- 
alized to the case where the right-hand side of equation ([6]) has the form 

kf kg 
fe=l k=l 

with fk and g\ as in Remark |2. II 

4 Regularity up to the parabolic boundary 

The following result gives conditions on the data which are sufficient for Holder continuity on 

[o,r] x n. 

Theorem 4.1 Let a G (0, 1), T > 0, JV > 2, and £1 C be a bounded domain with C 2 -smooth 
boundary Y. Let the assumptions (Hl)-(HS) be satisfied. Suppose further that 

u G B£"(p), g G Y D := B^^\[0, T}; L p (T)) n L p ([0, T]; Bpp * (r)) 
wt/j p > i + y, and t/iat i/ie compatibility condition 

uq = g\t=o on T 

is satisfied. Then for any bounded weak solution u of in VLt such that u — g a.e. on (0, T) x T, 
there holds 

McV-([o,nxU) - c (l u li=o(o T ) + |«o| + |/U r( [o,T];i 9 (n)) + |s|y D ) (13) 

-fpp (,_« »J 

urafft positive constants e = e(A, 2/, a,p, r, q, AT, f2) and C = C(A, i/, a,p, r, q, N, f2, T). 
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Proof. By Theorem 12. 31 the problem 

d?{v - u ) - Av = 0, t G (0,T), a; £ 
v = g, te(0,T),xGT, 
V\ t=0 = MO, X G n, 

admits a unique strong solution v in the class 

« G Z :=H"([0,T};L p (n))nL p ({0,T};H 2 p (n)) 

and 

Mz < Co(KI 2 _jl + |ff|y D ), 
s pp p ° (O) 

where Co only depends on a,p, N, T, f2. As in the proof of Theorem l3 . 1 1 we see that v G C" 5 ([0, T] x 
Q) for some S > 0. Furthermore, the mixed derivative theorem implies that 

A« 6 H"#([0,T];i p (O)) n L p ([0,T\;H$(n)) ^ fl^ ([0,T]; ff^(fi)) 

for all ? £ [0, 1]. Without restriction of generality we may assume that p £ (— + y, ~ + -^0- 
With 

ap p 

and <r := ^ - i G (0,1) we then have iZ^ ([0, T]; H^(n)) ^ L 2p (n T ), which shows that 

|A«| 2 G L p {Vl T ) with i + § < 1. 

Setting id = u — v, w is a bounded weak solution of 

d?w - DiiaijDjw) = / + I ),:<!,,!),<■> - Av, t G (0, T), at G ft, 

and w = a.e. on (0, T) x I\ 

Next, let f2o be an arbitrary bounded domain containing Q such that dist(fi, dClo) > 0. 
We extend w, /, ay and := DjV to [0, T] x fi by setting w,f,ipi — and Oy = Sij on 
[0, T] x (fi \ ft). Then w solves 

<9 t a u> - Di(aijDjw) = f + Difajipj - ipt), t e (0,T), x G Oo, 

in the weak sense, and thus Theorem 13. II and Remark 13.11 imply that w is Holder continuous on 
[0, T] x ft. Since u = v + w, the assertion of Theorem 14.11 follows . □ 

5 Proof of the global solvability theorem 

The proof of Theorem II . H is divided into three parts, devoted respectively to local well-posedness, 
existence of a maximally defined solution, and to a priori estimates which lead to global existence. 
Recall that the data belong to the following regularity classes: 



/ G X T := L p {[0,T];L p {ft)), u G Y 1 := B pp *"* {ft) 
geYS :^B pp [1 ^^\[0,T};L p (T))nL p ([0,T};B 2 p ;-(T)). 
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We seek a unique solution u of (JXJ) in the space 

Z T :=i/;([0,T];L p (fl))nl p ([0,T];ff p 2 (O)). 

1. Local well-posedness. Short-time existence and uniqueness in the regularity class Z s 
can be established by means of the contraction mapping principle and maximal L p -regularity for 
an appropriate linearized problem. We proceed similarly as in [27], see also [2] and |20) . 

We first define a reference function w G Z T as the unique solution of the linear problem 

d"(w - u ) - a i j(uQ)D i D j w = / + ay(«o)A«o DjU , t G (0,T), x e fl, 

w = g, te(0,T),xGT, 
w\ t =o =u , ie a, 

see Theorem 12.31 Note that the condition p > N + — ensures the embedding 
and thus we also have 

Z T C([0,T];F 7 ) cao.TjjC 1 ^)). 
For 5 G (0, T] and p > let 

£(£, p) = {w G Z s : v\ t=0 = uo, \v - w\ z s < p}, 

which is a closed subset of Z s . By Theorem 12.31 we may define the mapping F : £(<$, p) — > Z s 
which assigns to u E £(<5, p) the unique solution v = F(u) of the linear problem 

d?(v - u ) - a,ij(u )DiDjV = / + h(u, Du, D 2 u), t G (0, 6), x G Q, 

i' = g, t £ (0,6), xeT, (14) 
v|t=o = «0j a: 6 ^, 

where 

h(u, Du, D 2 u) — (a,ij(u) — aij(uo))DiDjU + a'tj(u)DiU DjU. 

Observe that every fixed point u of F is a local solution of ([1} and vice versa, at least for some 
small time interval [0,6]. 

Since Z s ^ C([0, 5]; C 1 ^)) we may set 

fJ-w(S) :— max{|iy(t, x) — uq(x)\ + \Dw(t, x) — Duq(x)\ : t G [0, 6], x G 

Evidently, p w (S) — > as 5 — > 0, due to iw|{=o = u . Letting u G X(<5, /?) we then have for all 
t G [0, 6} and cc G U 

\u(t,x) - u (a;)| + \Du(t,x) - Dw (x)| < |tt- w\ c([0 S] . cl(U)) +fi w (6) 

< M \u-w\ z s + p w (S) < M p + fi w (6), (15) 
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where the embedding constant Mo > does not depend on u and 8 £ (0,T]; the latter is true 
since u — w belongs to the space oZ 5 := {<p £ Z s : <p\t=o = 0}. (p~5|) yields for any u £ Y,(8,p) 
the bound 

\u{t,x)-u [x)\ + \Du{t,x)-Du (x)\<M po + fJ.w{T), t £ [0, 8], x £ U, (16) 

where we assume p £ (0,po]. 

Let now u £ £(£, p) and v = F(u). Then v — w £ qZ & solves the problem 

d"(v — w) — aij(uo)DiDj(v — w) = h(u, Du, D 2 u) — a' i; j(uo)DiUo DjUo, t £ (0, 5), x £ fl, 

v - w = 0, t £ (0,6), x £ r, 
(v — w)|f = o = 0, x £ il. 

Consequently, it follows from Theorem 12. 31 that for some constant Mi > which is independent 
of S £ (0, T] there holds 

\ v ~ w \z s ^ Mi\h(u, Du, D 2 u) — a' i j(uo)D i UQ DjUo\ x s 

< M-i\{aij(u) - aij(uo))DiDju\ X 5 + M^a'^ujD^ D 3 u - o^(u )A«o-Dj«o|x«- 
Using (fT5|) and (fT6|) we may estimate the first term as follows. 

\(aij(u) - aij(uo))DiDju\xs < (\A(u) - A(w)\^ L ^ N 2 + \A(w) - A(u )\^ L ^ )N 2) 

x (\D 2 u - D 2 w\ {x5)N 2 + \D 2 w\ (xS)N 2) 
< M 2 (p + p w (5))(p+\D 2 

W \(X s ) n2 J' 

where M 2 > does not depend on 6 and p. Similarly we obtain 

la'ij^DiuDjU - ay(u )Auo-Dj <uo\x* < M 3 {p + p w {S)) (p + 8?), 

with M3 > being independent of S and p; here the factor Sp comes from the estimate \z\ x s < 
(Ajv(0)5) l^l,*,. We conclude that 

\v-w\ z s <M((p + p(6)f, (17) 

where M and p(8) are constants, which do not depend on p, M is independent of 8, and p(8) is 
non-decreasing with p(S) — > as 8 — > 0. 

Next let Ui £ £(<5, p) and Vi = F(ui), i — 1,2. Then v\ — v 2 £ o% S solves the problem 

df(vi — v 2 ) — ciij(uo)DiDj(v — w) = h(ui, Dui, D 2 ui) — h(u 2 , Du 2 , D 2 u 2 ), t £ (0, 8), x £ Q, 

vi - v 2 = 0, t £ (0, 6), x £ r, 
(Vl - W2)|t=o =0, x e O, 

hence 

- «2|z^ < Mi|/i(ui, Dui, L> 2 Mi) - h(u 2 ,Du 2 , D 2 u 2 )\ x &. 
Estimating similarly as above we obtain 

\01-v2\z* < M((p + p(6))\ Ul -wlz', (18) 
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where M and fi(8) are like those in (|T7|) . 

Finally, the estimates (fT7|) and (fl8l) show that for sufficiently small p and 8 the mapping F 
is a strict contraction which leaves the set p) invariant. Local existence and uniqueness of 
strong solutions to ([1} follows now by the contraction mapping principle. 

2. The maximally defined solution. The local solution u G obtained in the first part 
can be continued to some larger interval [0,6 + Si] C [0,T]. In fact, let ug := u\ t= s G T 7 and 
define the set 

£(<J, 5i,p) := {w G 2" 5+Al : u| [0 ,5] = «, |« - < p}, 

where the reference function w G ^ T is now defined as the solution of the linear problem 

9" (to - uo) - a ij (u s )D i D J w = f + hi + X(s,T](t)a' ij (us)D i u s Djus, t G (0, T), a; G fi, 

to = g, t G (0,T), x G T, 
to| t =o = wo, ^ G fi, 

with 

^i = X[o,i5] (*)((<%'(«) - (Hj(u5))DiDjU + a' i:j (u)DiU D 3 uj. 

Observe that w|[o,a] = u, by uniqueness. So £(<5, <5i,p) is not empty and it becomes a complete 
metric space when endowed with the metric induced by the norm of Z s+Sl . 

Define next the mapping F : £(£, Si,p) — > Z 5+Sl which assigns to u G £(<5, 6i, p) the solution 
t> = -F(m) of the linear problem 

0f (w - tt ) - dijiu^DiDjV = f +~h(u,Du,D 2 u), t G (0, 5 + 5i), x G ft, 

w = g, te (0,6 + Si), x e r, 

to|t=o = t*0) a; G ft, 

where 

h(u,Du,D 2 u) = (a,ij(u) — ctij(us))DiDjU + a[j(u)DiuDjU. 

Since fi|[o,*l = u we have also t> | [o,<5] = u , by uniqueness. 
Observe that hi = h(u, Du, D 2 u) on [0,8] and thus 

[v-w\ z s+s 1 < Mi[h(u,Du,D 2 u) - a' tj (us)D t us D J us[ Lp{ [ S ^+s 1 }y.n)- 

Further, 

\F(ui) - F(u2)\ z s+H < Mi\h(ui,Dui,D 2 ui) - h(u 2 , Du 2 , D 2 u 2 )\L p ([s,s+8 1 ]xn), 

for Ui,tt2 G Ti(5,8i, p). Therefore we may estimate analogously to the first step to see that for 
sufficiently small 8i and p we have F(E(8, Si, p)) C £(<5, Si, p) and F is a strict contraction. Hence 
the contraction mapping principle yields existence of a unique fixed point of F in £(<5, Si, p), which 
is the unique solution of (JTJ) on [0, 8 + 8i]. 

Repeating this argument we obtain a maximal interval of existence [0, T max ) with T max < T , 
that is T max is the supremum of all r G (0, T) such that the problem (fT]) has a unique solution 
U G Z r . 

3. A priori bounds and global well-posedness. In order to establish global existence 
we will show that |u|z t stays bounded as t / T max . 
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Let r G (0,T max ) and u G Z T be the unique solution of (Q]). Setting bij(t,x) = a,i 3 ;(u(t, x)), it 
is evident that w is a weak solution of 

d?(u - uo) - I ),[b, , I = /, t G (0, r), x e ft. 

Since Ky C(ft) and Yg <^-» C([0,r] x T), Theorem 12.11 implies a uniform sup-bound for u|, 
namely 

\u(t,x)\ < d, t G [0,t], a; G ft, 

where the constant C*i depends only on the data |/|x T > Iffloo, woloo, 7\ o;, JV, and not on r. 
It follows then from Theorem 14.11 that for some e > we have 

\ U \c c ([0,r]xn) - C*2, 

where the number C% > 1 depends only on |/|x T i Ifflyji I^oIk,? T, a, N, and z/, not on r. In 
particular, we obtain a uniform Holder estimate for the coefficients bij, i, j = 1, . . . , N. 
The first equation of (|TJ) can be rewritten as 

d"(u — Uo) — bijDiDjU = f + a' i: j(u)DiuDjU. 

By Theorem 12. 31 the linear problem 

d?(v - u ) - hjDiDjV = f, te (0, r), x G ft, 
v = g, te (o,t), ier, 
v\ t =o = u , x e ft, 

has a unique solution v G i? T and there exists a constant Mi > independent of r such that 

\u - v\ Z r < M^a'^ujDiU Dju\ x ^ 

N 

<M X £ max l^-^lll^l 2 !^. (19) 

The assumption on p implies p > y and thus 

^(ft)^#i p (ft)^C £ °(ft) 

for some e € (0, e]- By the Gagliardo-Nirenberg inequality, there exists then 9 G (0, |) such 
that 

l-*M*> OI^^R-v, < CT|«(t, -)l!r»(n)l«(*i 01% < C^K*, -)|^ (n)) * G [0,r], 
and hence by Holder's and Young's inequality 

||-Dm| 2 | xt < C\Du\ 2 L2p{[Q ^ T]xn . RN) < C 3 \u\f p ( [0 iT ] ; h2(q))T p 

<c 4 H|t <£iH^ + c 5 (£i,^C4), 

for all £i > 0. This together with (fT9]) yields a bound for |u — which is uniform w.r.t. r. 
Since stays bounded as r / T max , it follows that |u|z t enjoys the same property. Hence 

we have global existence. □ 



14 



6 Decay estimate 

In this section we prove an Li decay estimate for solutions of ([TJ with / — and g = 0, that is 
we consider 

d?(u - m ) - Di{a,ij{u)Dju) = 0, t > 0, x e CI, 

u = o,t>o,ier, (20) 

u\ t =o = wo, x e O. 



Theorem 6.1 Let 57 C R N (N > 2) be a bounded domain with C 2 -smooth boundary T. Let 

2- — 

a G (0, 1), p > N + — , and uq € -Bpp p ° (f2) smc/i that uq — on T. Assume that condition 
(Q2) is satisfied. Then for the global strong solution u of i20\) the function {t i-> \u(t, Olz, 3 (o)l 
is continuous on [0, oo) and we have 

\,.(t \\* < C| " ol ^(") *>n 

wt/j positive constants c — c(a) and fi — n(y, N, O). 

Proof. Let T > be arbitrary and it G Z T be the solution of (l2Tfl) on [0, T]. We multiply the first 
equation in (|20)) by u, integrate over f2, and integrate by parts. Using the Dirichlet boundary 
condition this yields 

:= / ( udfu + a,ij {u)DjuDiU J dx = / uuogi- a if) dx =: I^if), a.ai£(0,T). 
Jn ^ 'hi 

Thanks to (Q2) and Theorem E3 with H = L 2 (ty we have for a.a. t £ (0,T) 

/i(t) > 5^1^,01^(0)+ 5 Si-a(*)|«(*, 011,(0) +H^(* ) -)li a( n ; R*)- 
On the other hand, Young's inequality implies that 

h{t)< \gi- a {t)\u{t,-)\l 2(n) + iffi-a(t)|uoli a( n)> a.a.te(0,T). 
Combining these estimates gives 

OlLfO) + MDu{t, Oli 2( o ;R «) < <?i-«(*)Kl! 2 (o), a.a. * e (0,T), 
which in turn, by Poincare's inequality, implies 

d?Ht,-)\l 2{n) +ti\u(t,-)\l 2 ( n) < 9i- a (t)\u \l 2in) , a.aJe(0,T), (21) 
where /i = p(v, N, 51) is a positive constant. Setting 

W(t) = \u(t, 0|| 2( o) and W = W(0) = | Uo |! 2 (o) 
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the fractional differential inequality (|21l) is equivalent to 

d?(W-W ) + fj,W <0 a.e. on(0,T). (22) 
Next, let V denote the solution of the corresponding fractional differential equation, that is 
d?(V -V )+ (iV = a.e. on(0,T), V(0) = V o = W o . (23) 
By the comparison principle for linear fractional differential equations (cf. [5]), we then have 

W(t)<V(t), te[0,T}. 

The solution of (f2"3"]l is given by 

V(t) = V E a (-fj,t a ), t£[0,T], 
where E a denotes the Mittag-LefHer function defined by 

fc=0 v ' 

see [TTJ Section 4.1]. Note that E 1 (z) = e z . It is known that for a G (0, 1) E a is a completely 
monotonic function in (— oo, 0] (see e.g. [TB]) and that there exists a constant c > such that 

E a (-x) < x>0, 
1 + x 

see [13j Formula (13)]. It follows that 

Ht,-)\l 2{n) < \u a \l 2m E a (-f,t a ) < 1 + L ^ ] , t G [0,T]. 
Since T > was arbitrary, the theorem is proved. □ 
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